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摘要 :【 目 的 】 蜜 蜂 是 典型 的 具有 发 育 狭 温 性 的 全 变态 昆虫 。 本 研究 以 对 低温 最 敏感 的 意大利 蜜蜂 
Apis mellifera ligustica 预 肾 为 研究 对 象 ,通过 低温 胁迫 不 同时 间 的 差异 表达 基因 (人 differentially 
expressed genes, DEGs) 趋势 分 析 , 探 讨 低温 胁迫 对 蜜蜂 发 育 影 响 的 关键 基因 。【 方 法 】 对 3 日 龄 意 
KFR AH E CFOuESSERT 207€ 低温 胁迫 18 h(T18) 和 36 h( T36) ,以 未 经 低温 胁迫 的 预 肾 为 对 照 
(CK) ,通过 Illumina HiSeq 平台 进行 转录 组 学 测定 。 利 用 Short Time-series Expression Miner 
(STEM) AF 3t 2 个 处 理 组 与 对 照 组 比较 共有 的 差异 表达 基因 进行 趋势 分 析 , 再 进一步 对 显著 富 
集 趋 势 模式 中 富 集 的 差异 表达 基因 进行 CO 分 类 和 KEGG pathway 分 析 。 利 用 RT-qPCR 对 随机 挑 
选 的 5 个 DEGs 的 表达 模式 进行 验证 。【 结果】 对 检测 到 1 062 个 T18 vs CK 和 T36 vs CK 共有 的 
DEGs 进行 趋势 分 析 , 发 现 3 个 显著 的 基因 表达 模式 ,包括 2 个 上 调 表达 模式 (Profile 6,74 539 个 基 
因 ;Profile 7, 有 271 个 基因 ) ,1 个 下 调 表达 模式 (Profile 1, 有 183 个 基因 ) 。 对 3 个 显著 富 集 趋势 模 
A DEGs 分 别 进行 GO 富 集 分 析 和 KEGG pathway 分 析 , 在 Profile 6 中 找到 同时 富 集 在 FoxO 信号 通 
路 和 寿命 调节 信号 通路 上 的 胰岛 素 样 肽 基因 ILP 和 又 头 蛋白 0 基因 FoxO 持续 上 调 ,说 明 低 温 胁 
迫 会 影响 蜜蜂 预 肾 晓 皮 激素 信号 传递 ;在 Profile 7 中 CREB 结合 蛋白 基因 CBP 持续 上 调 , 蜜 蜂 预 肾 
受到 低温 胁迫 会 影响 细胞 发 育 ; 在 Profile 1 中 细胞 色素 P450 基因 CYP450 306a1 (phm) 和 WNTI 显 
著 下 调 , 说 明 化 肾 时 受到 低温 胁迫 会 影响 晓 皮 激素 的 合成 。 通 过 RT-qPCR 分 析 挑 选 的 5 个 基因 的 
表达 结果 和 高 通 量 测序 结果 一 致 。【 结论 ] 通 过 对 蜜蜂 响应 低温 胁迫 的 差异 表达 基因 趋势 分 析 发 
现 ,蜜蜂 胰岛 素 和 晓 皮 激素 共同 调控 的 FoxO 可 能 是 低温 胁迫 抑制 蜜蜂 化 晴 的 一 个 关键 基因 。 本 
研究 为 探索 低温 胁迫 影响 蜜蜂 发 育 变 态 的 分 子 调 控 机 制 黄 定 了 基础 。 
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Trend analysis of differentially expressed genes in prepupae of Apis 


mellifera ligustica under low temperature stress 

YAO Dan, ZHOU Shu-Jing, ZHOU Bing-Feng, XU Xin-Jian' , ZHU Xiang-Jie" ( College of Animal 
Sciences ( College of Bee Science) , Fujian Agriculture and Forestry University, Fuzhou 350002, China) 
Abstract: [ Aims] Honey bees are of the holometabolous insects, whose brood is typically stenothermic. 
In order to understand the key genes involved in the sensitivity of the Italian honey bee, Apis mellifera 
ligustica to cold stress, we conducted a trend analysis of differentially expressed genes ( DEGs) of cold- 
sensitive prepupae in response to cold stress. [ Methods] The 3 d-old post-capped prepupae were exposed 
to cold temperature (20*C) for 18 (T18) and 36 h ( T36) , respectively, with those unexposed to cold 


stress used as the control ( CK), and then subjected to transcriptomic sequencing based on Ilumina 
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HiSeq™. Short Time-series Expression Miner ( STEM) software was applied to compare the shared DEGs 
among the three datasets. In addition, GO and KEGG pathway enrichment analyses were performed for 
DEGs with up- and down-regulation profiles. The expression profiles of five randomly selected DEGs were 
verified by RT-qPCR. [Results] Trend analysis showed that 1 062 DEGs shared by T18 vs CK and T36 
vs CK were clustered into three significant profiles including two up-regulation profiles ( Profile 6 covering 
539 DEGs and Profile 7 covering 271 DEGs) and one down-regulation profile ( Profile 1 covering 183 
DEGs). GO and KEGG pathway enrichment analyses showed that insulin-like peptide gene (ILP) and 
forkhead box protein O gene ( FoxO) were up-regulated in both FoxO signaling and longevity regulating 
pathways within Profile 6, suggesting that ecdysone signaling responds to cold stress during honey bee 
pupation. The CREB-binding protein gene ( CBP) in Profile 7 was up-regulated, and CYP450 306al 
(phm) and WNT| in Profile 1 were down-regulated, suggesting that these genes regulate cell 
development and ecdysone biosynthesis in response to low temperature. Besides, the expression profiles 
of the five selected genes based on RT-qPCR analysis were consistent with those based on the 
transcriptome. [Conclusion] The trend analysis of DEGs reveals that FoxO regulated by insulin signaling 
and ecdysone signaling in honey bees may act as a key gene that represses honey bee pupation in response 
to cold stress. Our data lay a foundation for exploring the molecular regulation mechanisms underlying 


metamorphosis in response to cold stress in honey bees. 


Key words: Apis mellifera ligustica; pupation; low temperature stress; transcriptomics; ecdysone 





蜜蜂 不 仅 为 人 类 提供 了 丰富 多 样 的 蜂 产品 ,而 
且 是 农业 中 最 重要 的 授粉 昆虫 ,在 维护 生态 平衡 发 
挥 重 要 作用 (Allen-Wardell et al., 1998) 。 近 年 来 ， 
全 球 气候 变化 和 病虫害 危害 导致 蜜蜂 的 种 群 数量 大 
幅 下 降 ( Milbrath et al., 2015) ,蜜蜂 蜂 群 的 减少 引起 
了 人 们 对 蜜蜂 健康 的 极 大 关注 。 

环境 温度 是 影响 昆虫 分 布 形态. 生理、 行为 等 
的 主要 生态 因素 之 一 (Burnett，1949; Ponsonby and 
Copland, 1996) 。 独 居 性 昆虫 可 通过 休 眼 或 滞 育 等 方 
式 耐 受 低温 等 不 良 环境 ,生态 幅 一 般 较 宽 (Tauber and 






































2017; Wang et al., 2017 ; Xu et al., 2017 ; Qin et al., 
2019), Æ OC F 20 H WS rp fE 5E ME Apis cerana 
cerana 成 蜂 转 录 组 测序 检测 到 与 环境 进程 相关 的 
Ca“ 信 和 号 通路 .Fox0 信号 通路 、Hippo 信号 通路 等 ， 
同时 热 激 蛋白 、 锌 指 蛋白 和 丝氨酸 / 苏 氮 酸 蛋白 激酶 
等 差异 表达 (Xu et al., 2017) 。 

相对 其 他 日 龄 的 蜜蜂 封 盖 子 ,3 日 龄 封 盖子 处 
于 幼虫 - 肾 变 态 的 过 渡 期 ,表现 出 对 低温 胁迫 最 敏 
感 。 前 期 研究 发 现 ,不 同日 龄 封 盖子 经 20Y 低温 胁 
迫 后 3 日 龄 封 盖子 最 早出 现 大 量 封 盖子 死亡 ,羽化 













































































Tauber, 1976) 。 而 蜜蜂 的 卵 . 虫 . 晴 ( 统 称 为 蜂 子 ) 是 
在 蜂 群 维持 较 稳 定 的 集 温 下 发 育 的 ,因此 形成 了 狭 温 
性 特点 ,发 育 温 区 为 29 ~ 38% ,最 适 发 育 温度 为 
35*C (Tautz et al., 2003; Groh et al., 2004; 朱 翔 杰 
SE, 2006) 。 偶 离 最 适 温度 ,羽化 后 成 蜂 的 外 部 形态 
(Medrzycki et al., 2010) 、 翅 脉 发 育 ( 周 冰峰 等 ,2011 ; 
Zhu et al., 2018) JHAR fk (Groh et al., 2004) 、 
寿命 (Wang et al., 2016) 和 行为 (Jones et al., 2005; 
Becher et al., 2009) 等 均 受 到 影响 。 蜜 蜂 发 育 对 温 
度 的 敏感 性 ,使 得 蜂 子 成 为 研究 温度 对 个 体 发 育 影 
响 的 理想 材料 ,有 助 于 认识 温度 与 生命 的 关系 。 
低温 对 昆虫 的 影响 主要 表现 在 糖 代 谢 、 脂 类 代 
谢 、 氮 基 酸 代谢 、 味 叭 代谢 和 硫 胶 素 代谢 等 代谢 通路 
上 富 集 大 量 差异 表达 基因 (differentially expressed 
genes，DEGs)( 库 尔 班 . 吐 松 等 ,2016; Tang et al., 




































































蜜蜂 的 寿命 也 最 短 ( Wang et al., 2016) 。 到 目前 为 
止 , 尚 无 研究 对 蜜蜂 封 盖子 发 育 期 间 低温 敏感 的 机 
制 给 出 解释 。 为 了 深入 研究 低温 对 蜜蜂 化 晴 影 响 的 
分 子 机 制 ,本 研究 采用 转录 组 学 方法 对 低温 处 理 不 
同时 间 的 意大利 蜜蜂 Apis mellifera ligustica 3 日 龄 封 
盖子 的 共有 差异 表达 基因 进行 趋势 分 析 , 从 分 子 水 
平 探讨 蜜蜂 化 肾 时 受 低温 影响 的 关键 基因 和 关键 通 
路 。 本 研究 对 丰富 昆虫 生态 学 和 昆虫 发 育 生物 学 有 
要 意义 。 


















































1 材料 与 方法 
1.1 供 试 蜜蜂 

样本 取 自 春季 增长 阶段 的 意大利 蜜蜂 A m. 
ligustica 蜂 群 ,通过 限制 蜂王 产 卵 空间 ,获得 卵 龄 一 
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Si B St EDD, 3 d 后 将 卵 脾 放 入 同一 个 哺育 群 哺 
育 。 移 入 哺育 群 5 d 后 这 些 幼 虫 即将 封 盖 , 获 取 4h 
内 同时 封 盖 的 封 盖子 ( 封 盖 梨 房 内 的 幼虫 和 晴 ) , 逢 
取样 本 梨 脾 , 放 入 恒温 恒 湿 箱 (35 + 0. 2 ，RH 
75% ) | CTHI-250B， 施 都 凯 仪 器 设备 (上 海 ) 有 限 公 
司 ] 中 培育 至 3 日 龄 (对 照 组 ,CK)。 再 将 一 部 分 3 
日 龄 封 盖 子 放 入 温度 20 上 0.2%C RH 75% ( Wang et 
al., 2016) 的 恒温 恒 湿 箱 中 分 别处 理 18 h 和 36 h 
(分 别 标记 为 T18 和 T36) ,作为 低温 处 理 组 。 所 取 
的 封 盖子 样本 立即 用 液 氮 冷冻 , 放 入 - 80C 保存 备 
用 。 每 个 样本 用 30 只 封 盖 子 (3 个 蜂 群 ,10 头 / 群 ) 
混合 建 库 和 高 通 量 测序 。 
1.2 高 通 量 测序 及 测序 数据 质量 评估 

利用 Trizol 法 提取 1.1 节 样 品 全 虫 总 RNA, 用 
人 带 有 Oligo (dT) W R R E R mRNA 后 ,加 入 
fragmentation buffer 使 其 片段 成 为 短片 段 的 mRNA 
为 模板 。 用 六 碱 基 随 机 引物 ( random hexamers ) 合 
成 cDNA 第 1 链 , 并 加 入 缓冲 液 dNTPs, RNase H 和 
DNA polymerase I 合成 cDNA 第 2 链 ,经 过 QiaQuick 
PCR 试剂 盒 纯 化 并 加 EB. 缓冲 液 洗 脱 经 末端 修复 、 
加 碱 基 A ,加 测序 接头 ,再 经 琼脂 糖 凝 胶 电泳 回收 目 
的 大 小 片段 ,并 进行 PCR 扩 增 ,从 而 完成 整个 文库 
制备 工作 ,构建 好 cDNA. 文库 ,高 通 量 测序 委托 广州 
基 迪 奥 生物 科技 有 限 公 司 完成 ,测序 平台 为 umina 
HiSeq™。 对 下 机 数据 的 raw reads 进行 过 滤 ,去 除 含 
adapter 的 reads, N 比例 大 于 1096 的 reads 和 低 质 
量 的 reads , 最终 得 到 clean reads ,然后 利用 Tophat2 
(2.1.1) 软件 将 reads 比 对 到 西方 蜜蜂 Apis mellifera 
的 参考 基因 组 (ncbi_GCF_000002195.4) 上 ,并 利用 
Cufflinks 组 装 转 录 本 ,转录 组 原始 数据 已 上 传 NCBI 的 
SRA 数据 库 (获取 号 : SRR10912937 - SRR10912945) 。 
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1.3 差异 表达 基因 趋势 分 析 

使 用 DESeq2 软件 ( V1. 20.0) 对 有 生物 学 重复 
的 处 理 组 与 对 照 组 间 (T18 vs CK, T36 vs CK) 基 因 
表达 量 进 行 差异 分 析 (Love et al., 2014) , DEGs ffi 
选 条 件 为 FDR( 错误 发 现 率 ) «0. 05，|log, Fold 
changel (Fold change, FC, 指 基因 表达 水 平 的 变化 
倍数 , FC = TACK, 其 中 了 为 处 理 组 的 基因 表达 水 
F, CK 为 对 照 组 的 基因 表达 水 平 ) 三 1。 利 用 STEM 
(Short Time-series Expression Miner) ,选择 参数 | -pro 
20 -ratio 1.0[ log,(2) =1, log; (1. 5) 20.5849625, 
log, (1. 2) = 0. 2630344 ] | ( Ernst and Bar-Joseph , 
2006) ,对 2 个 比较 组 间 共 有 的 差异 表达 基因 进行 趋 
势 分 析 , 在 软件 中 设置 P<0.05 为 差异 显著 的 表达 
模式 。 
1.4 ”差异 表达 基因 GO 及 KEGG pathway 富 集 分 析 

对 各 个 趋势 中 的 基因 进行 CO( Gene Ontology ) 
和 KEGG ( Kyoto Encyclopedia of Genes and Cenomes ) 
代谢 通路 进行 功能 注释 (Omicshare tools, http: // 
www. omicshare. com/tools/home/index/index. html ) , 
并 通过 假设 检验 计算 得 到 P 值 。 得 到 的 P 值 通 过 
FDR 校正 之 后 ,以 P 值 <0.05 为 国 值 ,满足 此 条 件 
的 GO 条目 (term) 和 通路 (pathway ) 定义 为 在 该 趋 
势 中 显著 富 集 的 GO term 和 Pathway。 
1.5 RT-qPCR 基因 验证 

随机 挑选 5 个 趋势 显著 的 基因 进行 RT-qPCR 
验证 , b-actin 为 内 参 基因 (引物 见 表 1)。 利 用 RNA 
抽 提 试剂 盒 提取 样本 总 RNA ,利用 反 转 录 试 剂 盒 合 
成 cDNA 第 1 链 (北京 全 式 金 生 物 技 术 有 限 公 司 ) 。 
以 上 述 的 cDNA 作为 模板 进行 qPCR, 反 应 体系 (20 
uL): H,O 6.6 uL, 2 x TransStart SuperMix 10 pL, 
上 下 游 引物 (5 umol/L)4& 1 uL, cDNA 模板 1 ul, 








































































































X1 RT-qPCR 引物 信息 
Table 1 Primers for RT-qPCR 

















基因 CenBank 登录 号 描述 引物 序列 (5 -3^) 
Cene CenBank accession number Gene description Primer sequence 
Actin 相关 蛋白 11 F: TGCCAACACTCTCCTTTCTG 
p-actin 413144 ; ; 
Actin-related protein 11 R: AGAATTGACCCACCAATCCA 
Kibra 和 蛋白 F: CGCCTATCAGCGAAACACCT 
kibra 413440 
Protein kibra R: GGCCTCGAATACTCCCGAAT 
j 40 F: TACCAGTGCTCCAGATCGAACG 
FoxO 727091 XXEH j 
Forkhead box protein O R: CGGCGCTACTGTTGCTTTCT' 
aip "DES 类 组 蛋白 乙酰 转移 酶 p300 F: TGTCCTGTATGCCTATCTTT 
f Histone acetyltransferase p300-like R: CGCTCTTCTCATTTCGGTTT 
新 肌 球 蛋白 IXb F: GAGGCACATCTTCCTTAGCG 
Myo20 551706 HIREA 
: Unconventional myosin-IXb R: CGTCACTCCTATCGTCGTCTT 
Armadillo 体 节 极 性 蛋白 F: ACCAGGGTGACCTGAGTTTC 
arm 408399 
Armadillo segment polarity protein R: GCCTCACGGATTCGTTTCTC 
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Dye 110.4 uL, qPCR 反应 条 件 : 95C 预 变性 4 min; 
95*C ZEE 15 s, 60% 退火 和 延伸 45 s, 共 40 个 循环 ; 
熔 解 曲线 默认 系统 程序 (QuantStudio5 ,美国 ) 。 样 
本 取 自 1.1 节 , 对 照 组 和 处 理 组 的 封 盖子 均 来 自 3 
个 蜂 群 ,每 个 蜂 群 取 1 头 进 行 混 合 提取 总 RNA。 每 
个 样本 设置 3 次 目的 基因 的 重复 扩 增 。 
1.6 数据 分 析 

目的 基因 的 RT-qPCR 表达 水 平 按照 内 参 基 因 
的 2 2 法 进行 计算 , 即 ACE = Ctagas 一 Cityssa 
( Livak and Schmittgen, 2001) ,各 基因 的 相对 表达 量 
经 log,(X) 转 换 后 与 高 通 量 测序 中 的 表达 水 平 进行 
比较 。 应 用 SPSS V21.0 软件 对 每 个 基因 人 处理 组 与 

















对 照 组 之 间 平 均 相 对 表达 量 的 差异 显著 性 进行 非 配 
对 了 检验 。 
2 结果 


2.1 测序 结果 质量 
通过 对 CK, T18 fü T36 测序 ,分别 得 到 原始 读 
段 平均 为 56 759 342 , 45 977 301 和 54 196 127 ,过 
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Profile 6: 539 genes 





Profile 7: 271 genes Profile: 1:183 genes Profile 0: 64 genes 


EA 


滤 后 得 到 高 质量 有 效 读 段 平均 分 别 为 55 578 624, 
45 032 615 和 53 121 097。 统 计 两 端的 Q20 和 Q30 
分 别 平均 在 98. 20% 和 94. 85% 以 上 上。 表明 本 研究 
测序 数据 质量 良好 ,可 用 于 进一步 分 析 。 
2.2 差异 表达 基因 趋势 

通过 差异 表达 基因 (DECs ) 分 析 ,Tl18 vs CK 有 
1 394 个 DEGs,T36 vs CK 有 2 722 个 DEGs ,说 明 随 
着 低温 处 理 的 时 间 增 长 , 受 影响 的 基因 数 明 显 增加 。 
其 中 两 比较 组 共有 DEGs 1 062 个 ,Tl8 vs CK 特有 
DEGs 有 332 个 ,T36 vs CK 特有 DEGs 有 1 660 个 ， 
上 调 趋势 的 DEGs 显著 多 于 下 调 趋 势 的 DECs。 

对 两 比较 组 共有 的 1 062 个 DEGs 进行 趋势 分 
析 , 在 8 个 基因 表达 模式 中 有 3 个 为 显著 基因 表达 
模式 (P<0.05), 包 括 2 个 上 调 表 达 模 式 , Profile 6 
模式 为 基因 表达 量 先 增加 后 保持 不 变 ,有 539 个 基 
;Profile 7 模式 为 基因 表达 量 随 低温 胁迫 时 间 延 
长 一 直上 调 , 有 271 个 基因 ;1 个 下 调 表达 模式 
Profile 1 ,为 随 低温 胁迫 时 间 延 长 基因 表达 量 先 下 调 
表达 ,之 后 保持 不 变 , 有 183 个 基因 (图 1)。 多 数 基 
因 归 入 上 述 3 个 显著 基因 表达 模式 。 
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图 1 低温 胁迫 不 同时 间 3 日 龄 
Fig. 1 





意大利 蜜蜂 封 盖子 差异 表达 基因 表达 趋势 
Trend analysis of differentially | genes ( DEGs) in the 3 d-old post-capped brood of 





Apis mellifera ligustica subjected to low temperature stress for different time 
A; DEGs 的 趋势 分 析 结果 ,趋势 类 型 方 框 内 数字 为 己 值 The result of trend analysis of DEGs, the number in the box denoting the P-value; B; Profile 
6 包含 DEGs 的 表达 模式 Expression profile of DEGs in Profile 6; C; Profile 7 包含 DEGs 的 表达 模式 Expression profile of DEGs in Profile 7; D; 








Profile 1 包含 


DEGs 的 表达 模式 Expression profile of DEGs in Profile 1. CK: 对 照 组 (未 经 低温 胁迫 的 3 日 龄 封 盖 子 ) Control group (3 d-old post- 


capped brood unexposed to low temperature) ; T18 : 20% 低温 处 理 18 h 的 3 日 龄 封 盖子 3 d-old post-capped brood exposed to low temperature (20*C ) 
for 18 h; T36: 20% 低温 处 理 36ph 的 3 日 龄 封 盖 子 3 d-old post-capped brood exposed to low temperature (20% ) for 36 h. 下 图 同 The same below. 


2.3 差异 表达 基因 的 GO 分 类 

GO 富 集 分 析 结 果 显 示 ,3 日 龄 意 蜂 封 盖子 在 受 
到 低温 胁迫 后 DEGs 归 类 到 生物 学 进程 的 基因 数 最 
多 。Profile 6 趋势 上 的 DEGs 分 别 富 集 在 30 个 GO 
Z HE, Profile 7 趋势 上 的 DEGs 富 集 在 28 个 GO 








条 目 上 。 这 两 种 趋势 从 生物 学 进程 分 析 , 主要 富 集 
在 细胞 进程 .代谢 进程 . 单 有 机 体 进 程 ;从 分 子 功能 
分 析 , 主要 富 集 在 结合 .催化 活性 、 分 子 传感器 活性 

从 细胞 组 件 分 析 , 主 要 富 集 在 细胞 膜 组 分 、 细 胞 膜 、 
细胞 。 除 此 之 外 Profile 6 趋势 中 的 DEGs 还 富 集 在 
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核酸 结合 转录 因子 活性 生物 学 进程 负 调控 发 育 进 
FE ` 突 触 和 突 触 组 分 上 ,而 Profile 7 趋势 中 还 富 集 在 
行为 (1 个 )、 转 录 因 子 活 性 、 蛋 白质 结合 (1 个) M 
胞 连接 等 条 目 上 (图 2) 。 

Profile 1 为 随 低温 胁迫 时 间 延 长 基因 表达 量 先 
下 调 表达 , 之 后 保持 不 变 ,该 趋势 上 富 集 183 个 
DEGs ,这 些 基因 寅 集 在 25 个 GO 条 目 上 , 排 在 前 10 
的 GO 条 目 分 别 为 代谢 进程 (28 个 ) 、 单 有 机 体 进程 


细胞 成 分 Cellular component 
膜 封 闭 腔 Membrane-enclosed lumen 


(26 个 ) .结合 (22 个 ) ` 细 胞 进程 (21 个 )、 催 化 活性 
(18 个 ) 细胞 膜 组 分 (12 个 ) ` 细 胞 膜 (12 个 ) ` 细 胞 
(12 个 ) .细胞 组 分 (12 个 ) 应 激 (8 个 ) 条 目 上 (图 
2) 。 说 明 低温 胁迫 主要 影响 3 日 龄 封 盖子 的 代谢 
进程 和 细胞 结构 等 。 同 时 机 体 核酸 结合 转录 因子 活 
性 ,发 育 相 关 的 基因 不 会 随 着 低温 胁迫 时 间 延 长 表 
达 量 改变 ,与 行为 相关 的 基因 表达 量 随 着 低温 胁迫 
时 间 延 长 而 升 高 。 
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—— Profile 7 
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转录 因子 活性 、 蛋 白质 结合 Transcription factor activity, protein binding 
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结构 分 子 活 性 Structural molecule activity 
分 子 功能 调节 器 Molecular function regulator 


转运 子 活性 Transporter activity 


信号 转运 子 活性 Signal transducer activity 

核酸 结合 转录 因子 活性 Nucleic acid binding transcription factor activity 
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2 低温 胁迫 不 同时 间 3 日 龄 意大利 蜜蜂 封 盖 子 差异 表达 基因 GO 富 集 分 析 
Fig. 2 GO enrichment analysis of differentially expressed genes in the 3 d-old post-capped brood 





of Apis mellifera ligustica subjected to low temperature stress for different time 
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2.4 差异 表达 基因 的 KEGG pathway 富 集 分 析 
通过 KEGG pathway 富 集 分 析 结 果 显 示 , Profile 
6 趋势 上 的 DEGs 富 集 在 70 条 信号 通路 上 ,其 中 有 
48 条 与 代谢 相关 的 信号 通路 。 基 因数 富 集 最 多 的 
为 代谢 途径 (22 个 ) 、 神 经 活性 受 体 (11 个 ) KER 
谢 物 的 生物 合成 (11 个 ) 不同 环境 中 的 微生物 代谢 
(8 个 ) .抗生素 的 生物 合成 (6 个 ) 、 细 胞 外 基质 - 受 
体 相 互 作用 (5 个 ) 氨基 酸 的 生物 合成 (5 个 ) DERE 
体 (5 个 )、 胞 知 作 用 (5 个 ) 等 ;找到 基因 ILP 
(insulin-like peptide receptor gene, ID411297 ) 和 基 
FoxO( forkhead box protein O gene, ID 727091 ) 同 
时 富 集 在 FoxO 信号 通路 和 寿命 调节 信号 通路 上 
( 表 2)。 

Profile 7 趋势 中 的 271 个 DEGs 富 集 在 51 条 信 
号 通路 上 ,其 中 与 代谢 相关 的 信号 通路 有 31 条 。 基 
因数 富 集 最 多 的 为 代谢 进程 (16 个 ) 次 生 代谢 物 的 
生物 合成 (7 个 )、 味 叭 代谢 (4 个 ) 、 匮 氨 酸 降解 (3 












































































































































个 ) Jak-STAT 信号 通路 (2 A) Notch 信号 通路 (2 
个 ) 淀粉 和 蔗糖 代谢 (2 个 ) TGF-B 信号 通路 (2 
个 ) 基础 转录 因子 (2 个 ) Hippo 信号 通路 (2 个 ); 
其 中 E CBP ( CREB-binding protein gene, 
ID726332) 同时 富 集 在 Jak-STAT, Notch, TGF-p, 
FoxO 和 Wnt 5 条 信号 通路 上 ( 表 3)。 

Profile 1 模式 中 183 个 DEGs 富 集 在 43 条 信号 
通路 上 。DEGs 主要 富 集 在 代谢 进程 模块 中 ,包括 
新 陈 代谢 总 览 中 的 代谢 途径 (13 个 ) 不 同 环境 中 微 
生物 代谢 (3 个 ) 次 生 代 谢 物 的 生物 合成 (3 个 ) ; 核 
ERR PA RIR G 个 ) 和 喀 喧 代谢 (3 个 ); 
氨基 酸 代谢 中 有 酷 损 酸 代谢 (2 个 ) WAAR RILA 
酸 和 和 谷 氮 酸 代谢 (2 个 ) ;有 1 个 基因 CYP450 306a1 
(phm, cytochrome P450 306al gene，ID08398 ) 富 集 
在 昆虫 激素 的 生物 合成 信号 通路 上 ,还 发 现 基因 
WNTI(protein Wnt-1 gene, ID 2j 413502) 同时 富 集 在 
Hippo, Wnt 和 hedgehog( Hh)3 条 信和 号 通路 上 ( 表 4) 。 







































































x 




















表 2 Profile 6 中 DEGs 基因 数 富 集 前 12 的 KEGG 通路 
Table2 Top 12 KEGG pathways enriched by DEGs of Profile 6 
















































































通路 基因 数 通路 ID P fÑ 

Pathway Number of genes Pathway ID P-value 
代谢 途径 Metabolic pathway 22 ko01100 0.9426 
神经 活性 受 体 Neuroactive receptor 11 ko04080 0. 0000 
次 生 代 谢 物 的 生物 合成 Biosynthesis of secondary metabolites 11 ko01110 0.3507 
不 同 环境 中 的 微生物 代谢 Microbial metabolism in diverse environments 8 ko01120 0. 1227 
抗生素 的 生物 合成 Biosynthesis of antibiotics 6 ko01130 0. 6883 
细胞 外 基质 - 受 体 相互 作用 Extracellular matrix -receptor interaction 5 ko04512 0.0019 
氨基 酸 的 生物 合成 Biosynthesis of amino acids 5 ko01230 0.0553 
拼接 体 Spliceosome 5 ko03040 0.3436 
METEJ Endocytosis 5 ko04144 0.3784 
Hippo 信号 通路 Hippo signaling pathway 4 ko04391 0. 1053 
FoxO 信和 号 通路 FoxO signaling pathway 3 ko04068 0.2902 
寿命 调节 信和 号 通路 Longevity regulating pathway 2 ko04213 0.4714 

表 3 Profile7 中 DEGs 基因 数 富 集 前 12 的 KEGG 通路 
Table3 Top 12 KEGG pathways enriched by DEGs in Profile 7 

通路 基因 数 通路 ID P (H 

Pathway Number of genes Pathway ID P-value 
代谢 进程 Metabolic process 16 ko01100 0.3555 
次 生 代 谢 物 的 生物 合成 Biosynthesis of secondary metabolites 7 ko01110 0. 2070 
味 叭 代谢 Purine metabolism 4 ko00230 0. 1691 
赖 氨 酸 降解 Lysine degradation 3 ko00310 0.0357 
Jak-STAT 信号 通路 Jak-STAT signaling pathway 2 ko04630 0. 0466 
Notch 信和 号 通路 Notch signaling pathway 2 ko04330 0.0554 
淀粉 和 蔗糖 代谢 Starch and sucrose metabolism 2 ko00500 0.0902 
TGF-B 信号 通路 TGF-beta signaling pathway 2 ko04350 0.1124 
基础 转录 因子 Basal transcription factors 2 ko03022 0.1239 
Hippo 信和 号 通路 Hippo signaling pathway 2 ko04391 0.2370 
FoxO 信和 号 通路 FoxO signaling pathway 1 ko04068 0.6254 
Wnt 信和 号 通路 Wnt signaling pathway 1 ko04310 0.6874 
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5&4 Profile 1 rh DEGs 基因 数 富 集 前 13 的 KEGG 通路 
Table4 Top 13 KEGG pathways enriched by DEGs in Profile 1 

通路 基因 数 通路 ID P fH 
Pathway Number of genes Pathway ID P-value 
代谢 途径 Metabolic pathway 13 ko01100 0.0500 
WE EFC Pyrimidine metabolism 3 ko00240 0.0462 
喇叭 代谢 Purine metabolism 3 ko00230 0.1337 
内 质 网 的 蛋白 质 加 工 Protein processing in endoplasmic reticulum 3 ko04141 0. 1384 
不 同 环境 中 微生物 代谢 Microbial metabolism in diverse environments 3 ko01120 0.1891 
次 生 代 谢 物 的 生物 合成 Biosynthesis of secondary metabolites 3 ko01110 0.5612 
酷 氨 酸 代谢 Tyrosine metabolism 2 ko00350 0.0127 
丙 氨 酸 、 天 冬 氨 酸 和 谷 氨 酸 代 谢 Alanine, aspartate and glutamate metabolism 2 ko00250 0.0375 
Hippo 信号 通路 Hippo signaling pathway 2 ko04391 0.1028 
核糖 体 Ribosome 2 ko03010 0.3623 
昆虫 激素 的 生物 合成 Insect hormone biosynthesis 1 ko00981 0.1175 
Wnt 信号 通路 Wnt signaling pathway 1 ko04310 0.4957 
Hh 信和 号 通路 Hedgehog signaling pathway 1 ko04340 0.2217 





2.5 差异 表达 基因 的 RT-qPCR 验证 
利用 RT-qPCR 技术 检测 随机 挑选 的 5 个 DEGs 
(kibra, FoxO, CBP, Myo20 和 arm) 的 表达 模式 。 
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结果 显示 这 些 DEGs 的 表达 水 平 与 RNA-Sed 测序 数 
据 的 基因 表达 水 平 变化 趋势 一 致 (图 3 ) ,说 明 本 次 
转录 组 的 测序 数据 真实 可 靠 。 
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图 3 低温 胁迫 不 同时 间 3 日 龄 意大利 蜜蜂 封 盖子 差异 表达 基因 RNA-Seq 测序 数据 的 RT-qPCR 验证 





Fig. 3 Verification of RNA-Seq sequencing data of differentially expressed genes in the 3 d-old post-capped brood 


of Apis mellifera ligustica subjected to low temperature stress for different time by RT-qPCR 
A: T18 vs CK; B: T36 vs CK. 图 中 数据 为 平均 值 上 标准 差 , 柱 上 双星 号 代表 同一 基因 处 理 组 与 对 照 组 之 间 相 对 表达 量 差异 极 显 著 ( 忆 <0.01， 


非 配对 了 检验 ) Data in the figure are mean € SD, and the double asterisk above bars represents extremely significant difference between the treatment 














group and the control group of the same gene in the relative expression level ( P 0.01, unpaired T-test). 


3 讨论 


真 社 会 性 昆虫 蜜蜂 蜂 子 具有 典型 的 狭 温 性 发 育 
特点 , 蜂 群 维持 稳定 的 合适 集 温 对 蜂 子 的 健康 发 育 
具有 十 分 重要 意义 。3 日 龄 封 盖子 是 蜜蜂 进行 化 肾 
起 始 阶段 ( 预 师 期 ) ,其 低温 敏感 性 强 。 本 研究 采用 
趋势 分 析 对 两 个 梯度 时 间 的 低温 胁迫 预 师 转 录 本 表 
达 模 式 进 行 聚 类 分 析 ,得 到 了 梯度 低温 胁迫 条 件 下 























的 DEGs 及 归 类 ,筛选 出 具有 生物 学 意义 的 基因 集 
和 蜜蜂 响应 低温 胁迫 的 关键 基因 。 在 此 基础 上 ,我 
们 重点 讨论 了 FoxO 信和 号 通路 中 的 fox0 上 调 表 达 
在 低温 胁迫 后 蜜蜂 预 晴 无 法 晓 皮 化 师 中 的 可 能 
机 制 。 

低温 胁迫 可 能 导致 旷 皮 激素 (20- 
hydroxyecdysone, 20E) 的 合成 减少 。20E Æ Bj E% 
化 合 物 ,是 在 昆虫 前 胸腺 (prothoracic glands, PG) 中 
以 胆固醇 或 植物 固 醇 类 物质 为 底 物 ,经 一 系列 酶 , 特 
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别 是 矶 22-5 4L | x 2-8 15 RE, xe 25 AER ME 
化 作用 合成 的 (Gilbertet et al., 2003; Gilbertet, 
2004) 。 在 对 下 调 模式 Profile 1 中 的 DEGs 进行 
KEGG 富 集 分 析 显 示 ,在 激素 合成 信号 通路 上 存在 
一 个 下 调 基因 CYP450 306al (phm) , 它 通过 调节 碳 
25- 凑 化 酶 的 活性 来 参与 20E 合成 ,phm 表达 量 下 调 
将 会 导致 碳 25- 凑 化 酶 活性 降低 ,进而 影响 合成 20E 
(Niwa et al., 2004; Warren et al., 2004) 。 因 此 推测 
低温 会 影响 20 E 合成 ,导致 其 滴 度 降低 。 

低温 胁迫 后 ,3 日 龄 封 盖子 无 法 化 晴 , 可 能 是 由 
于 ILP 调控 FoxO 磷酸 化 留 在 细胞 质 内 ,无 法 行使 其 
功能 而 导致 无 法 晓 皮 。 昆 虫 化 晴 时 ILP 调控 FoxO 
完成 赔 皮 的 过 程 如 下 :在 昆虫 中 存在 状 椎 动物 的 胰 
岛 素 类 似 物 ILP, 它 的 重要 靶 标 蛋白 为 Fox0 , fe ILP 
上 调 时 ,Fox0 发 生 磷 酸化 ,无 法 进入 细胞 核 而 留 在 
细胞 质 内 ,促进 细胞 生长 , 延 开 化 师 (Pan et al., 
2018) ; 当 ILP 下 调 时 ,Fox0 去 磷酸 化 和 人 核 与 20E 和 
EcR-Usp 聚合 体 激 活 后 ,诱导 Br-C, E75A 和 E75B 
等 初级 及 DER3 , DHRA 和 fz- 有 1 等 次 级 应 答 基 因 表 
达 , 来 激活 凑 肽 酶 -A ( CPA ) 表达 完成 赔 皮 发 生 
(Conzy et al., 2002; Zhou and Riddiford, 2002; Wu 
and Brown, 2006; Minakuchi et al., 2008; Jinder et 
al., 2013)。 在 研究 中 上 调 模 式 Profile 6 的 差异 表 
达 基 因 的 富 集 分 析 结 果 显 示 , 在 晓 皮 相关 的 Fox0 
信号 通路 (3 个 ) 与 寿命 相关 的 寿命 调节 信和 号 通路 
(2 个 ) 同 时 富 集 到 了 ILP 和 FoxO 基因 。ILP 的 表 
达 量 上 调 会 导致 Fox0 磷酸 化 后 无 法 入 核 留 在 细胞 
质 内 ,无 法 被 20E 和 EcR-Usp 聚合 体 激活 ,导致 20E 
的 初级 应 答 基因 Br-C 无 法 接受 20E 信号 , 旷 皮 信和 号 
传递 受阻 。 另 外 ,20E 具有 持 抗 胰岛 素 的 类 似 物 ILP 
的 作用 ,同时 与 ILP 共同 调控 FoxO 信号 通路 ,从 前 
一 段 论述 低温 胁迫 可 能 导致 20E 滴 度 降低 ,也 正好 
与 该 结果 相符 。 因 此 推测 在 3 日 龄 封 盖子 在 受到 低 
温 胁迫 后 , Fox0 发 生 磷 酸化 留 在 细胞 质 内 不 能 入 
核 ,是 阻碍 蜜蜂 化 师 的 一 个 关键 基因 。 化 师 是 一 个 
极其 复杂 的 过 程 , 暗 皮 只 是 其 中 的 一 个 环节 ,除了 蜡 
皮 受 到 影响 之 外 ,低温 胁迫 带 来 的 影响 应 该 还 有 很 

我 们 前 期 研究 结果 显示 3 日 龄 意 蜂 封 盖子 在 受 
到 低温 胁迫 后 寿命 显著 缩短 (Wang et al., 2016), 
FoxO 在 果 晶 的 寿命 调节 起 到 重要 作用 (Tatar et al., 
2004) 。 转 录 组 数据 显示 FoxO 基因 在 受到 低温 胁 
迫 时 表达 量 发 生 改 变 ,这 也 可 能 是 受到 低温 胁迫 后 ， 
蜜蜂 寿命 缩短 的 原因 之 一 。 


























































































































































































































致谢 福建 农林 大 学 科研 助手 刘 也 齐 参与 样本 获取 
工作 ,硕士 研究 生 刘 一 名 、 冯 害 营 和 村寨 参与 部 分 文 
献 研 究 工 作 ,本 科 生 冯 海 成 、 藉 心 怡 、 赫 旦 畅 ' 常 斐然 
参与 部 分 实验 ,在 此 表示 感谢 ! 


参考 文献 (References) 


Allen-Wardell G, Bernhardt P, Bitner R, Burquez A, Buchmann S, 
Cane J, Cox PA, Dalton V, Feinsinger P, Ingram M, Inouye D, 
Jones CE, Kennedy K, Kevan P, Koopowitz H, Medellin R, 
Medellin-Morales S, Nabhan GP, Pavlik B, Tepedino V, Torchio 
P, Walker S, 1998. The potential consequences of pollinator 
declines on the conservation of biodiversity and stability of food crop 
yields. Conserv. Biol., 12(1) : 8 - 17. 

Becher MA, Scharpenberg H, Moritz RFA, 2009. Pupal developmental 
temperature and behavioral specialization of honeybee workers ( Apis 
mellifera L. ). J. Comp. Physiol. A, 195(7) : 673 - 679. 

Burnett T, 1949. The effect of temperature on an insect host-parasite 
population. Ecology, 30(2) : 113 — 134. 

Ernst J, Bar-Joseph Z, 2006. STEM a tool for the analysis of short time 
series gene expression data. BMC Bioinform., 7(1) 191. 

Gilbert LI, 2004. Halloween genes encode P450 enzymes that mediate 
steroid hormone biosynthesis in Drosophila melanogaster. Mol. Cell. 
Endocrinol., 215(1 —2) ; 1 - 10. 

Gilbert LI, Rybczynski R, Warren JT, 2003. Control and biochemical 
nature of the ecdysteroidogenic pathway. Annu. Rev. Entomol., 47; 
883 - 916. 

Gonzy G, Pokholkova GV, Peronnet F, Mugat B, Demakova OV, 
Kotlikova IV, Lepesant JA, Zhimulev IF, 2002. Isolation and 
characterization of novel mutations of the Broad-Complex, a key 
regulatory gene of ecdysone induction in Drosophila melanogaster. 
Insect Biochem. Mol. Biol., 32(2) : 121 - 132. 

Groh C, Jürgen Tautz, Wolfgang Róssler, 2004. Synaptic organization in 
the adult honey bee brain is influenced by brood-temperature control 
during pupal development. Proc. Natl. Acad. Sci. USA, 101 
(12) : 4268 - 4273. 

Jinder M, Palli SR, Riddiford LM, 2013. The juvenile hormone 
signaling pathway in insect development. Annu. Rev. Entomol., 58; 
181 - 204. 

Jones JC, Helliwell P, Beekman M, Maleszka R, Oldroyd BP, 2005. 
The effects of rearing temperature on developmental stability and 
learning and memory in the honey bee, Apis mellifera. J. Comp. 
Physiol. A, 191(12) : 1121 - 1129. 

Livak KJ, Schmittgen TD, 2001. Analysis of relative gene expression 


2-7^59 method. 


data using real-time quantitative. PCR and the 
Methods , 25(4) : 402 - 408. 

Love MI, Huber W, Anders S, 2014. Moderated estimation of fold 
change and dispersion for RNA-seq data with DESeq2. Genome 
Biol., 15(12) ; 550. 

Medrzycki P, Sgolastra F, Bortolotti L, Bogo G, Tosi S, Padovani E, 


Porrini C, Sabatini A, 2010. Influence of brood rearing temperature 


on honey bee development and susceptibility to poisoning by 


696 昆虫 学 报 Acta. Entomologica Sinica 63 卷 





pesticides. J. Apic. Res., 49(1) : 52 — 59. 

Milbrath MO, van Tran T, Huang WF, Solter LF, Tarpy DR, Lawrence 
FK, Huang ZY, 2015. Comparative virulence and competition 
between Nosema apis and Nosema ceranae in honey bees ( Apis 
mellifera). J. Invert. Pathol., 125; 9 - 15. 

Minakuchi C, Zhou X, Riddiford LM, 2008. Krüppel homolog 1 ( Kr- 
hl) mediates juvenile hormone action during metamorphosis of 
Drosophila melanogaster. Mech. Dev., 125(1 -2) : 91 —105. 

Niwa R, Matsuda T, Yoshiyama T, Namiki T, Mita K, Fujimoto Y, 
Kataoka H, 2004. CYP306A1, a cytochrome P450 enzyme, is 
essential for ecdysteroid biosynthesis in the prothoracic glands of 
Bombyx and Drosophila. J. Biol. Chem., 279 (34); 35942 - 
35949. 

Pan J, Di YQ, Li YB, Chen CH, Wang JX, Zhao XF, 2018. Insulin 
and 20-hydroxyecdysone oppose each other in the regulation of 
phosphoinositide-dependent kinase-l expression during insect 
pupation. J. Biol. Chem., 293 (48): 18613 - 18623. 

Ponsonby DJ, Copland MJW, 1996. Effect of temperature on 
development and immature survival in the scale insect predator, 
Chilocorus nigritus ( F. ) ( Coleoptera; Coccinellidae ). Biocontrol 
Sci. Technol., 6(1): 101 - 110. 

Qin M, Wang H, Liu Z, Wang Y, Zhang W, Xu B, 2019. Changes in 
cold tolerance during the overwintering period in Apis mellifera 
ligustica. J. Apic. Res., 58(5): 1-8. 

Tang B, Liu XJ, Shi ZK, Shen QD, Xu YX, Wang S, Zhang F, Wang 
SG, 2017. Transcriptome analysis and identification of induced 
genes in the response of Harmonia axyridis to cold hardiness. Comp. 
Biochem. Physiol. Part D Genomics Proteomics , 22; 78 —89. 

Tatar M, Hwangbo DS, Gersham B, Tu MP, Palmer M, 2004. 
Drosophila dFOXO controls lifespan and regulates insulin signaling 
in brain and fat body. Nature, 429(6991) ; 562 —566. 

Tauber MJ, Tauber CA, 1976. Insect seasonality; —diapause 
maintenance, termination, and postdiapause development. Annu. 
Rev. Entomol., 21; 81 —107. 

Tautz J, Maier S, Groh C, Róssler W, Brockmann A, 2003. Behavioral 
performance in adult honey bees is influenced by the temperature 
experienced during their pupal. Proc. Natl. Acad. Sci. USA, 100 
(12) : 7343 — 7347. 

Tusong K, Lu XY, Liu XN, Ma J, 2016. Transcriptome analysis of 
short-term low temperature stress response of Carapax trionycis from 


Junggar. Acta Entomol. Sin., 59(6): 581 — 591. [gZKBE - üt 


松 ， 陆 雪 莹 , 刘 小 宁 , 马 纪 , 2016. ENS ZI e FR A SEC 
迫 响应 的 转录 组 分 析 . 昆虫 学 报 , 59(6) : 581 591] 
Wang J, Gao G, Zhang R, 2017. Metabolism and cold tolerance of 











Chinese white pine beetle, Dendroctonus armandi, ( Coleoptera: 
Curculionidae: Scolytinae) during the overwintering period. Agric. 
For. Entomol., 19(1) : 10 -22. 

Wang Q, Xu X, Zhu X, Zhou S, Huang Z, Zhou B, 2016. Low- 
temperature stress during capped brood stage increases pupal 
mortality, misorientation and adult mortality in honey bees. PLoS 
ONE, 11(5) : e0154547. 

Warren JT, Petryk A, Marqués G, Parvy JP, Shinoda T, Itoyama K, 
Kobayashi J, Jarcho M, Li Y, O'Connor MB, Dauphin-Villemant 
C, Gilbert L, 2004. Phantom encodes the 25-hydroxylase of 
Drosophila melanogaster and Bombyx mori; a P450 enzyme critical 
in ecdysone biosynthesis. Insect Biochem. Mol. Biol., 34 (9): 
99] — 1010. 

Wu Q, Brown MR, 2006. Signaling and function of insulin-like peptides 
in insects. Annu. Rev. Entomol., 51; 1 —24. 

Xu K, Niu Q, Zhao H, Du Y, Jiang Y, 2017. Transcriptomic analysis 
to uncover genes affecting cold resistance in the Chinese honey bee 
(Apis cerana cerana). PLoS ONE, 12(6) : e0179922. 

Zhou BF, Zhu XJ, Li Y, 2011. New variation of hind wing veins of Apis 
cerana cerana caused by low temperature. Acta Ecol. Sin., 31(5) : 
1387 -1392. [ 周 冰峰 , RIR, EH , 2011. 低温 导致 中 华 蜜 
蜂 后 翅 翅 脉 的 新 变异 . 生态 学 报 , 31(5): 1387 - 1392] 

Zhou X, Riddiford LM, 2002. Broad specifies pupal development and 
mediates the * status quo’ action of juvenile hormone on the pupal- 
adult transformation in Drosophila and Manduca. Development , 129 
(9) : 2259 - 2269. 

Zhu X, Xu X, Zhou S, Wang Q, Chen L, Hao Z, Zhou B, 2018. Low 
temperature exposure (20°C ) during the sealed brood stage induces 
abnormal venation of honey bee wings. J. Apic. Res. 57(3): 
458 —465. 

Zhu XJ, Zhou BF, Luo QH, Luo F, Zhong GX, Shi ML, 2006. Effects 
of temperature on development of bee worker bees in sealing period. 
New Agric. Technol., (Special Issue) : 57 - 60. [ 朱 翔 杰 , 周 冰 
E, P HAE, 罗芳 , 钟 贵 喜 ,史明 亮 , 2006. 温度 对 蜜蜂 工蜂 封 
盖 期 发 育 的 影响 . 农业 新 技术 ，( 特刊 ) : 57 -60] 





























(责任 编辑 : 马 丽 萍 ) 


